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Abstract 



We present evidence for the first observation of electromagnetically bound tt im- 
pairs (7rK-atoms) with the DIRAC experiment at the CERN-PS. The 7riir-atoms 
are produced by the 24 GeV/c proton beam in a thin Pt-target and the 7r =l= and K T - 
mesons from the atom dissociation are analyzed in a two- arm magnetic spectrometer. 
The observed enhancement at low relative momentum corresponds to the production 
of 173 ± 54 7rK-atoms. The mean life of 7ri^-atoms is related to the s-wave ttK- 
scattering lengths, the measurement of which is the goal of the experiment. From 
these first data we derive a lower limit for the mean life of 0.8 fs at 90% confidence 
level. 

Key words: DIRAC experiment, exotic atoms, scattering length, 7rK-scattering, 
chiral perturbation 

PACS: 36.10.-k, 32.70.Cs, 25.80.Nv, 29.30.Aj 



1 Introduction 



The study of electromagnetically bound hadronic pairs is an excellent method 
to probe QCD at very low energy. Opposite charge long-lived hadrons such 
as pions and kaons can form hydrogen-like atoms bound by the Coulomb 
interaction. The strong interaction leads to a broadening of the atomic levels 
and dominates the lifetime of the atom. 

Pion-pion interaction at low energy, constrained by the approximate 2-fiavour 
SU(2) (u,d) chiral symmetry, is the simplest and well understood hadron- 
hadron process [T1I21I3] . The observation of 7r + 7r~-atoms (A 2n ) was reported in 
ref. |1] and a measurement of their mean life in ref. [5]. 

The low energy interaction between the pion and the heavier (strange) kaon is 
a tool to study the more general 3-fiavour SU(3) (u,d,s) structure of hadronic 
interaction, which is not accessible in 7T7r-interactions. A detailed study of nK- 
interaction provides insights into a potential flavour (u,d,s) dependence of the 
crucial order parameter or quark condensate in Chiral Perturbation Theory 
(ChPT) ®. 

A measurement of the 7ri^-atom (A^k) lifetime was proposed already in 1969 
[7] to determine the difference |ai/ 2 — a 3 / 2 | of the s-wave 7ri^-scattering lengths, 

* This publication is dedicated to the memory of Ludwig Tauscher. 
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where the indices 1/2 and 3/2 refer to the isospin of the Tri^-system. The 
7r K^-atom decays predominantly by strong interaction into the neutral me- 
son pair 7T°K (n°K ). The decay width of the vrA'-atom in the ground state 
is given by the relation 08]: 



T(A nK ) = — = - a 3 fx 2 p* \a 1/2 - a 3/2 \ 2 (1 + 5). (1) 

tis y 

Tis is the lifetime of the atom in the ground state, a the fine structure constant, 
\x the reduced if^K^ mass, and p* = 11.8 MeV/c the outgoing K° or it 3- 
momentum in the nK center-of-mass system. The term 5 ~ (4 ± 2) % |8] 
accounts for corrections, due to isospin breaking and the quark mass difference 
fn u — rfid- Hence a measurement of T(A n K) provides a value for the scattering 
length \ai/2 — a.3/2 1 - The mean life of vri^-atoms is predicted to be 3.7 ± 0.4 fs 

The width T(A wK ) can also be determined from the s-wave phase shifts ob- 
tained from 7r.fr-scattering, i.e. from the interaction of kaons with nucleons. 
The s-wave phase shifts are, however, poorly known due to the absence of 
data below 600 MeV/c and, correspondingly, the uncertainties in ai/ 2 and a 3 / 2 
are substantial. The overall interaction is attractive (attractive in the isospin 
1/2 and repulsive in the 3/2 state). The Roy-Steiner equations lead, with the 
available scattering data, to results [H] that are neither consistent with the 
most precise measurements 1 1 Qj nor with predictions from ChPT |llj . 

The method used by DIRAC is to produce pions and kaons with a high en- 
ergy proton beam impinging on a thin target. Pairs of oppositely charged 
mesons may interact and form electromagnetically bound systems. Their sub- 
sequent ionization in the production target leads to mesons emerging from 
the target with low relative momentum (thereafter called atomic pairs). The 
mean life of the 7r-fT-atom can then be calculated from the number of observed 
low-momentum pairs. This method was first proposed in 1985 [12] and was 
successfully applied to 7r + 7r~-atoms in Serpukhov on the U-70 synchrotron in- 
ternal proton beam [3], and with DIRAC-I at the CERN-PS on beam line T8. 
Data from DIRAC-I lead to a mean life r 27r = (2.91±g g2) fs for 7r + 7T~-atoms 

® 

Cross sections for A^x-production have been calculated in ref. [13]. In this 
paper we report on the observation of 7ri^-atoms from the first data with 
DIRAC-II. 
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2 Experimental setup 



Details on the initial apparatus (DIRAC-I) used to study 7r + 7r~-atoms can 
be found in ref. [14]. A sketch of the modified spectrometer (DIRAC-II) used 
to collect the nK (and more tttt) data is shown in fig. [TJ The 24 GeV/c pro- 
ton beam (1) from the CERN-PS impinges on a 26 fim Pt-target (2). The 
spill duration is 450 ms with an average intensity of 1.6 x 10 11 protons/spill. 
The proton beam then passes through a vacuum pipe and is absorbed by 
the beam dump. The secondary particles are collimated through two steel 
shielding blocks (3) and (7), upstream of the microdrift chambers (4) and 
downstream of the ionization hodoscope (6), respectively. They pass through 
a vacuum chamber (8) and are bent by the 1.65 T field of the dipole magnet 
(9). The two-arm spectrometer is tilted upwards with respect to the proton 
beam by an angle of 5.7°. Positive particles are deflected into the left arm, 
negative ones into the right arm. 




Fig. 1. Sketch of the DIRAC-II spectrometer (top view). 1 - proton beam, 2 - pro- 
duction target, 3 - shielding, 4 - microdrift chambers, 5 - scintillation fiber detector, 
6 - ionization hodoscope, 7 - shielding, 8 - vacuum chamber, 9 - dipole magnet, 10 
- drift chambers, 11 - vertical hodoscope, 12 - horizontal hodoscope, 13 - aerogel 
Cerenkov modules, 14 - heavy gas Cerenkov detector, 15 - N2-Cerenkov detector, 
16 - preshower detector, 17 - absorber, 18 - muon scintillation hodoscope. The solid 
lines crossing the spectrometer arms correspond to typical ir^- and if T -trajectories 
from the ionization of 7rX-atoms in the production target. 

We now describe in more detail the upgrade which was performed to search 
for and study vri^-atoms [IS]. The upstream detectors (4-6), see fig. [TJ were 
either replaced or upgraded. However, they were not yet fully operational 
during data taking, and were therefore not used in the analysis presented 
here. The tracking is performed by drift chambers (10) which have a spatial 
resolution of 85 /zm. The vertical hodoscope (11) consisting of 20 scintillating 
slabs with a time resolution below 140 ps is used for timing. The horizontal 
hodoscope (12), made of 16 horizontal scintillating slabs, is used for triggering 
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by selecting oppositely charged particles with a vertical displacement smaller 
than 75 mm. 

The N 2 -Cerenkov detector (15) was already used previously to reject electrons 
and positrons. The refractive index is n = 1.00029 and the average number of 
photoelectrons N pe = 16 for particles with velocity [3 = 1. The inner part of 
the original container had to be cut to clear space for the two new Cerenkov 
detectors needed for kaon identification. Since the momenta of the two mesons 
originating from the breakup of the TrX-atoms are very small in the center-of- 
mass system, they have similar velocities in the laboratory system, and hence 
kaons are less deflected than pions. Typical trajectories are shown in fig. [T| 

The heavy gas C 4 F 10 -Cerenkov detectors in both arms (14) identify pions but 
do not respond to kaons nor (anti)protons [IB]. Four spherical and four flat 
mirrors each focus the light towards the phototubes. The alignment of the mir- 
rors was checked with a laser beam [IT]. To keep a constant refractive index 
of n = 1.0014, the gas has to be cleaned permanently with a complex recir- 
culating system [18j. The average number of photoelectrons is 28 for particles 
with (3 = 1. 

The aerogel Cerenkov detector (13) in the left arm identifies kaons and rejects 
protons [i~9f2"0] . Such a detector is required only in the left arm since the 
contamination from antiprotons in the right arm is small due to their low 
production rate. The detector consists of three modules. Two modules of 121 
each (refractive index n = 1.015) cover the relevant momentum range between 
4 and 8 GeV/c. The aerogel stacks are 42 cm high and are read out by two 
5"-Photonis XP4570 photomutipliers with UV windows. The typical number 
of photoelectrons is 10 for (3 = 1 particles. A third overlapping module with 
13£ aerogel and n = 1.008 covers the small angle region to reject protons with 
momenta above 5.3 GeV/c. 

The aerogel tiles are stacked pyramidally to increase the radiator thickness 
halfway between the photomutipliers and to compensate for light absorption. 
Due to the low light yield for the n = 1.008 module, and the strong UV- 
light absorption, we use a wavelength shifter. The aerogel tiles are sandwiched 
between Tetratex foils on which tetraphenylbutadiene (TPB) has been evap- 
orated [TH]. The typical number of photoelectrons is 4 - 8 for (3 = 1 particles. 

The preshower detector (16) provides additional electron/hadron separation 
in the offline analysis. A lead converter, typically 10 - 25 mm thick, is placed 
in front of a 10 mm thick scintillator. An additional converter/scintillator is 
installed to compensate for the drop of efficiency of the N 2 -Cerenkov detector 
at small angles, in the region covered by the aerogel and heavy gas detectors. 
The iron absorber (17) and the array of scintillation counters (18) are used 
to suppress muons. 
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Pairs such as e + e~, n + ix~ , ix~ K + or n + K~ are selected by a two-level trigger. 
For n~K + - and 7r + i^ _ -candidates the first level trigger requires no signal in 
the heavy gas detector of the K-avm and in the N 2 -Cerenkov detectors. The 
two tracks have to cross the same (or one of the two adjacent) slab(s) in 
each of the horizontal hodoscopes. The two trajectories being asymmetric (see 
fig. [TJ, the kaon is required to cross slabs of the vertical hodoscope located 
in front of the heavy gas detector, while the associated pion has to fly at 
large angles in the opposite arm. The second level trigger uses raw hits 
from the drift chambers and, based on a lookup table, rejects pairs with high 
relative momentum. The accepted trigger rate is limited by the buffer memory 
of around 2000 events per spill. 



3 Tracking and calibration 

As mentioned already, only detectors downstream of the dipole magnet are 
used here for event reconstruction. The trajectories are determined by the drift 
chambers, the pattern recognition starting from the horizontal x-coordinate in 
the last plane and extrapolating back to the target. A straight line is first fitted 
through the hits and extended into the magnet yoke. A deflection algorithm 
[22] calculates the slope and x-coordinate of the track at the magnet entrance, 
following the magnetic field map. The trajectory is then extrapolated linearly 
to the target, with the constraint that the track origin has to coincide with the 
center of the beam spot. This determines the momentum of the particle. For 
the vertical y-coordinate the straight line from the drift chamber information 
is extrapolated back to the center of the beam spot at the production target 

pan]. 

The variable of interest in the following analysis is the relative momentum Q 
of the i^ ± 7r :F -pairs in their center-of-mass systems. In the transverse plane, the 
resolution on the relative momentum Q T (typically 3 MeV/c) is dominated by 
multiple scattering, while the resolution on the longitudinal component Ql 
(< 1 MeV/c) is not affected [23125] . For further analysis we use therefore only 
Ql- 




1.112 1.113 1.114 1.115 1.116 1.117 1.118 

jc"p invariant mass (GeV/c 2 ) 

Fig. 2. 7r~p-mass distribution in the A-region. The line shows the Gaussian fit. 
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The momentum calibration was cross-checked with tracks from A — > 7i~p 
decays. Figure [2] shows the invariant 7r~p-mass distribution. A Gaussian fit is 
applied leading to a mass of 1115.35 ± 0.08 MeV/c 2 (statistical error) and a 
width of a = 0.58 ± 0.01 MeV/c 2 , dominated by momentum resolution. 

We have verified that the new hardware and software were able to reproduce 
the signal from 7r + 7r~-atoms. Details of the procedure for 7r + 7r~ atoms can be 
found in ref. [51123]. Figure [3] shows the momentum distribution after back- 
ground subtraction. The enhancement at low Ql corresponds to 7098 ± 533 
atomic 7r + 7r~-pairs. Bound 7r + 7r~-pairs are not expected above \Q L \ = 2 MeV/c 
and indeed the distribution in fig. [3] is flat and compatible with zero, which 
validates the background subtraction method. Note that the present data can- 
not be compared directly with those of ref. [5] because we used here only the 
downstream detectors, and a different event selection was performed. 




2 4 6 8 10 12 14 

IQJ (MeV/c) 



Fig. 3. Ql- distribution measured with DIRAC-II for part of the 7r + 7r~ data af- 
ter background subtraction. The accumulation of events at low Ql is due to 
7r + 7r~-atoms. The curve is a Gaussian fit to guide the eye. 



4 Data analysis 

Figure [4] shows the four mechanisms which contribute to the production of 
tt^K^ -pairs. Accidental pairs are due to particles produced on different nucle- 
ons (fig. [4^), non-Coulomb-pairs are associated with the production of long- 
lived intermediate states (fig. [4)b) . On the other hand, 7r ± A' =F -pairs which in- 
teract electromagnetically form correlated Coulomb-pairs (fig. ^Sp), or atomic 
bound states (fig. The N A atoms, while traveling through the target, can 
either decay, be (de)-excited or break up into n A 7r ± A' =F -pairs which emerge 
from the target with very low relative momentum. 

We now describe the analysis steps [23]- For prompt pairs the time difference 
between the positive and negative spectrometer arm lies between -0.5 and 0.5 
ns. Accidental pairs are first removed using the time information from the 
vertical hodoscopes. Accidental pairs are also needed for subsequent analysis 
and we select those pairs with a time difference between -12 and -6 ns. The 
choice of the negative sign avoids the contamination from slow protons. The 
events have then to satisfy the following criteria: 
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a) b) c) d) 



Fig. 4. Production mechanisms of irK-pairs: a) accidental-pairs from two protons; 
b) non-Coulomb-pairs from long-lived intermediate states such as the 77-meson; c) 
Coulomb-pairs from direct production or from short-lived intermediate states; 
d) irK-aioms. 

• no electrons nor positrons, 

• no muons, 

• one drift chamber track per arm, 

• \Q L \ < 20 MeV/c, 

• Q T < 8 MeV/c, 

• the momentum of the kaon lies between 4 GeV/c and 8 GeV/c, 

• the momentum of the pion lies between 1.2 GeV/c and 2.1 GeV/c. 

With the excellent time resolution of the vertical hodoscope pions, kaons and 
protons below 2.5 GeV/c can be separated by time-of- flight [26]. For the ir~K + 
analysis the aerogel detector is used in addition to remove protons in the 
positive arm, while for the n + K~ analysis, the time difference between the 
negative and the positive arm measured with the vertical hodoscope has to 
be negative in order to remove protons faking pions. Once the accidentals 
have been subtracted the prompt pairs (N pr ) are composed of the following 



three types: atomic-pairs 
(N nC ). Therefore 



n 



Coulomb-pairs (N c ), and non-Coulomb-pairs 



JSP" = n A + N c + N 



jnC 



(2) 
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Fig. 5. a) QL-distributions for the ir~ K + data sample (26 fim Pt-target). The 
histogram shows the data. The fitted non- Coulomb- and Coulomb-pairs are in green 
(cl) and blue (c2) respectively, together with the total background in turquoise 
(c3). b) Residuals between data and fitted background. The solid line illustrates the 
distribution of atomic-pairs (see text). 

However, the background arising from 7r + 7r~- and 7r~p-pairs with misidenti- 
fied particles must be considered since the kaon flux is much lower than the 
pion and proton fluxes. Pions can be selected with the heavy gas detector in 
coincidence and protons with the aerogel detector in anticoincidence. We then 
determine Ql by assigning to the pion (or the proton) the mass of a kaon. For 
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7r + 7r~ events this incorrect mass assignment shifts the (^-distribution by - 
150 MeV/c [23] and therefore does not overlap with the Qi-distribution of ttK 
events. The contribution from 7r~p-pairs, non-Coulomb 7r~lT + -pairs and acci- 
dentals have a similar linear (^-distribution. We assume that the background 
due to Coulomb uncorrelated pairs can be described by the (^-distribution of 
accidentals, following a similar analysis for 7r + 7r~-atoms [5]. The non-Coulomb 
background also includes the background from 7r~p-pairs. Coulomb correlated 
pairs have to be simulated [27112%] . 

To determine the contribution from Coulomb- and non-Coulomb pairs we se- 
lect the momentum range 3 < \Ql\ < 20 MeV/c, where no atoms are expected. 
The (^-distribution is 



race 



dN pr dN c dN° 

dQ L dQ L dQ L 

where dN c /dQ L and dN acc /dQL are the (normalized) differential probabilities 
for Coulomb correlated or uncorrelated pairs, respectively. The fit variable f3 
is the corresponding number of correlated pairs. We choose a bin size of 0.25 
MeV/c. The x 2 -function to be minimized is 



X 



80 

E 

|i|=13 



dNP r 
dQL,, 



dN acc 

dQh.i 



a 



pr 



(4) 



where af r are the corresponding statistical errors in the measured number 
of prompt pairs. Figure [5^i shows the MINUIT results for Coulomb- and non- 
Coulomb contributions to n~K + events. Since the shapes of both contributions 
are known, one can extrapolate into the \Ql\ < 3 MeV/c signal region. The 
difference (residuals) between the data and the sum of both contributions is 
plotted in fig. [5}d. Above \Ql\ = 3 MeV/c the residuals are consistent with 
zero, while the enhancement at low relative momentum is the first evidence 
for 7r~i ; r + -atoms. 



A Gaussian distribution describes adequately the low momentum enhance- 
ment observed in the n + ir~ data [23]. To guide the eye we also apply a Gaus- 
sian fit here (line in fig. [5}d). The integral of the Gaussian distribution contains 
147 ± 61 atomic 7r~i^ + -pairs. 

A similar fit based on equ. §4§ is applied to ir + K~ events. However, the number 
of events is smaller, due to the lower production cross section for negative 
kaons. The fit results are summarized in table Q] The errors on the number of 
Coulomb-pairs (3 are the full (MINOS) errors, while the errors on n A are given 
by the square roots of the measured bin contents, statistical fluctuations on 
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±59 ±31 




1.24 


1356 


164 


29 


74 39 






±396 


±108 


±15 


±35 ±19 



Table 1 

Left: number f3 of Coulomb-pairs outside the signal region (3 < \Ql\ < 20 MeV/c), 
number N c of Coulomb-pairs extrapolated into the signal region (\Ql\ < 3 MeV/c), 
and number n A of detected atomic-pairs from the residuals of the fit. Right: expected 
number of atoms N A , calculated from the number of detected Coulomb-pairs, and 
expected number of atomic pairs n A using a breakup probability of 53%. 

the much larger Monte-Carlo sample being negligible. 

Figure [6] shows the sum of the n~K + and n + K~ residuals. We obtain 

n A {i^K^) = 173 ± 54 (5) 

detected atomic pairs with a statistical significance of 3.2a. The systematic 
uncertainty is estimated to be around 5%, much smaller than the statistical 
one. 



100 

£ 50 

5 

| 
Ul 

-50 
-100 

2 4 6 8 10 12 14 16 18 20 

IQJ (MeV/c) 

Fig. 6. Residuals between data and the fitted background for tt~ K + and tt + K~ . A 
Gaussian fit has been applied (solid line) to illustrate the distribution of atomic-pairs. 

The evidence for the observation of vr-fT-atoms is strengthened by the obser- 
vation of correlated Coulomb-pairs which, a fortiori, implies that atoms have 
also been produced. This can be seen as follows, without involving simulation: 
non-Coulomb pairs have a similar Qi-distribution as accidentals. Hence di- 
viding the normalized distribution for prompt pairs by the one for accidentals 
one obtains the correlation function R describing Coulomb-pairs. The func- 
tion R, shown in fig. [7] for tc~K + as a function of \Ql\, is clearly increasing 
with decreasing momentum, proving that Coulomb-pairs have been observed. 
In the signal region (\Ql\ < 3 MeV/c) one obtains 858 ± 247 Coulomb-pairs 
from the data in fig. [7j without resorting to Monte- Carlo. The same procedure 
can be applied to tx + K~ events, leading to 313 ± 148 Coulomb-pairs [2"3"] . 
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Fig. 7. Correlation function R as a function of \Ql\ for 7r~K + -pairs. The devi- 
ation from the horizontal dotted line proves the existence of Coulomb correlated 
7r~K + -pairs. 

The ratio k of the number of produced atoms to the number of Coulomb-pairs 
with small relative momenta has been calculated: k = 0.62 [T2"ll2"9"] . However, 
one needs to take into account the acceptance of the apparatus and the cuts 
applied in the analysis. By Monte-Carlo simulation |30] we determine the ratio 
k e xp = 0.24 between the number of atoms produced within acceptance and the 
number of detected Coulomb-pairs below \Q L \ = 3 MeV/c (and below Q T = 
8 MeV/c). This then leads to the expected number N A of atoms, 204 ± 59 
for n + K~ } and 74 ± 35 for tt~K + (table [TJ right). The uncertainty on k exp is 
negligible. 



The breakup probability 



P 



n; 



br 



(6) 



relates the number of atoms to the number of atomic pairs. A calculation of 
the breakup probability as a function of mean life (fig. [8j has been performed 
using the Born approximation [15]. For the predicted mean life of 3.7 ± 0.4 
fs [8j the corresponding breakup probability Pb r is 53% (dotted line in fig. |8~J. 
From equ. (|6j we then find the number n A of expected pairs given in table]TJ 
right. These numbers are in good agreement with the number n A of observed 
atomic-pairs (fifth column in table [TJ. 




T(fs) 



Fig. 8. Breakup probability Pb r for the 26 /im Pt-target as a function of mean life 
of 7rK-atoms in the ls-state. The horizontal solid line is the measured breakup 
probability and the horizontal dashed line the 1.28<r lower bound corresponding to 
a lower limit of 0.8 fs for the mean life. The excluded area (90% confidence level) is 
shown in turquoise. The horizontal dotted line gives the theoretical prediction |8J. 

Conversely, we use the number of observed atomic-pairs n A from the x 2 - 
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minimization and the number N c of Coulomb-pairs below \Ql\ < 3 MeV/c 
(table [l| to calculate the breakup probability P br from equ. ^ with N A = 
Kx P ■ N^- The result for tt^K^ (P br = 64 ± 25 %) is shown by the horizontal 
solid line in fig. |8j This leads to a lower limit for the mean life of 7ri^-atoms 
of ns = 0.8 fs at a confidence level of 90%. This result can be translated into 
an upper limit |ai/ 2 — 03/2! < 0.58 m~ l at 90% confidence level, in agreement 
with predictions [8f9] . 



5 Conclusions 



We have presented the first evidence for the production of 7ri^-atoms by de- 
tecting 173 ±54 atomic-pairs. The evidence is strengthened by the observation 
of correlated irK (continuum) Coulomb-pairs from which the number of bound 
states (atoms) is predicted and found to be in agreement with observation. A 
lower limit on the mean life of 0.8 fs is established with a confidence level of 
90%. We note that the choice of Pt as production target was driven by the 
high breakup probability facilitating the observation of 7rK-atoms. Data are 
now being collected for a more accurate measurement of the lifetime with e.g. 
a 98 Ni-target, for which the breakup probability is lower (~35% accord- 
ing to ref. [E]) but still rapidly rising around the predicted mean life of 3.7 fs. 
The ultimate goal of the experiment is to measure the lifetime of 7ri^-atoms 
with a precision of about 20% [15j. 
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